#set pond parameters
	pond=array(dim=c(15,6))	
	pond[1, ]=c(name='Ox1',      area=1029059,  bottom=1.97, surf=0,min.surf=0,max.surf=5)
	pond[2, ]=c(name='Ox2',      area=776105,  bottom=1.31, surf=5,min.surf=0,max.surf=5)
	pond[3, ]=c(name='Ox3',      area=161244,  bottom=0.66, surf=0,min.surf=0,max.surf=5)
	pond[4, ]=c(name='TM1',      area=79807,  bottom=0,    surf=0,min.surf=0,max.surf=5)
	pond[5, ]=c(name='TM2',      area=80111,  bottom=0,    surf=0,min.surf=0,max.surf=5)
	pond[6, ]=c(name='TM3',      area=53560,  bottom=0,    surf=0,min.surf=0,max.surf=5)
	pond[7, ]=c(name='Pilot',    area=65212,  bottom=0,    surf=0,min.surf=0,max.surf=5) 
	pond[8, ]=c(name='Chlorine', area=4169,  bottom=1.25, surf=0,min.surf=0,max.surf=5)
	pond[9, ]=c(name='Allen',    area=471875,  bottom=0,    surf=0,min.surf=0,max.surf=5) 
	pond[10,]=c(name='Gearheart',area=322845,  bottom=-.625,surf=0,min.surf=0,max.surf=5)
	pond[11,]=c(name='Hauser',   area=425062,  bottom=-1.88,surf=0,min.surf=0,max.surf=5)
	pond[12,]=c(name='Brackish', area=569500,  bottom=1,    surf=0,min.surf=0,max.surf=5)
	pond[13,]=c(name='Marsh Box',area=1000,  bottom=0,    surf=0,min.surf=0,max.surf=5)
	pond[14,]=c(name='Pond Box', area=1000,  bottom=0,    surf=0,min.surf=0,max.surf=5)
	pond[15,]=c(name='Bay',      area=10^10, bottom=0,    surf=0,min.surf=0,max.surf=5)

#set box parameters
	box=array(dim=c(39,9))
	#type 1=weir, 2=pump, 3=pipe, 4=sluice/valve
	box[1 ,]=c(type=1,weir.elev=4.93,top.elev=7.95,width=10,roughness=1,pump=0,upstream=1,downstream=2,name='Ox1 to Ox2')
	box[2 ,]=c(type=1,weir.elev=5.38,top.elev=7.82,width=10,roughness=1,pump=0,upstream=1,downstream=2,name='Ox1 to Ox2')
	box[3 ,]=c(type=1,weir.elev=5.27,top.elev=7.65,width=10,roughness=1,pump=0,upstream=2,downstream=4,name='Ox2 to T1')
	box[4 ,]=c(type=1,weir.elev=5.22,top.elev=7.61,width=10,roughness=1,pump=0,upstream=2,downstream=4,name='Ox2 to T1')
	box[5 ,]=c(type=1,weir.elev=5.22,top.elev=7.02,width=10,roughness=1,pump=0,upstream=2,downstream=5,name='Ox2 to T2')
	box[6 ,]=c(type=1,weir.elev=6.74,top.elev=7.01,width=10,roughness=1,pump=0,upstream=2,downstream=5,name='Ox2 to T2')
	box[7 ,]=c(type=1,weir.elev=0,top.elev=7.04,width=10,roughness=1,pump=0,upstream=2,downstream=3,name='Ox2 to Ox3')
	box[8 ,]=c(type=1,weir.elev=0,top.elev=7.07,width=10,roughness=1,pump=0,upstream=2,downstream=3,name='Ox2 to Ox3')
	box[9 ,]=c(type=1,weir.elev=0,top.elev=2.8,width=10,roughness=1,pump=0,upstream=4,downstream=13,name='T1 to T-Box')
	box[10,]=c(type=1,weir.elev=2.68,top.elev=3.75,width=10,roughness=1,pump=0,upstream=5,downstream=13,name='T2 to T-Box')
	box[11,]=c(type=1,weir.elev=2.71,top.elev=3.65,width=10,roughness=1,pump=0,upstream=5,downstream=13,name='T2 to T-Box')
	box[12,]=c(type=1,weir.elev=2.70,top.elev=3.88,width=10,roughness=1,pump=0,upstream=5,downstream=13,name='T2 to T-Box')
	box[13,]=c(type=1,weir.elev=2.66,top.elev=3.94,width=10,roughness=1,pump=0,upstream=5,downstream=13,name='T2 to T-Box')
	box[14,]=c(type=1,weir.elev=2.67,top.elev=3.94,width=10,roughness=1,pump=0,upstream=6,downstream=13,name='T3 to T-Box')
	box[15,]=c(type=1,weir.elev=2.64,top.elev=4.16,width=10,roughness=1,pump=0,upstream=6,downstream=13,name='T3 to T-Box')
	box[16,]=c(type=1,weir.elev=2.69,top.elev=4.27,width=10,roughness=1,pump=0,upstream=6,downstream=13,name='T3 to T-Box')
	box[17,]=c(type=1,weir.elev=2.69,top.elev=4.46,width=10,roughness=1,pump=0,upstream=6,downstream=13,name='T3 to T-Box')
	box[18,]=c(type=1,weir.elev=5.25,top.elev=6.94,width=10,roughness=1,pump=0,upstream=3,downstream=6,name='Ox3 to T3')
	box[19,]=c(type=1,weir.elev=0,top.elev=7.74,width=10,roughness=1,pump=0,upstream=3,downstream=14,name='Ox3 to P-Box')
	box[20,]=c(type=1,weir.elev=0,top.elev=7.74,width=10,roughness=1,pump=0,upstream=8,downstream=9,name='Cl to Allen')
	box[21,]=c(type=1,weir.elev=0,top.elev=7.2,width=10,roughness=1,pump=0,upstream=8,downstream=9,name='Cl to Allen')
	box[22,]=c(type=1,weir.elev=0,top.elev=7.78,width=10,roughness=1,pump=0,upstream=8,downstream=9,name='Cl to Allen')
	box[23,]=c(type=1,weir.elev=0,top.elev=7.7,width=10,roughness=1,pump=0,upstream=8,downstream=9,name='Cl to Allen')
	box[24,]=c(type=1,weir.elev=0,top.elev=7.61,width=10,roughness=1,pump=0,upstream=8,downstream=9,name='Cl to Allen')
	box[25,]=c(type=1,weir.elev=0,top.elev=7.26,width=10,roughness=1,pump=0,upstream=8,downstream=9,name='Cl to Allen')
	box[26,]=c(type=1,weir.elev=0,top.elev=7.15,width=10,roughness=1,pump=0,upstream=9,downstream=10,name='Allen to Gearheart')
	box[27,]=c(type=1,weir.elev=0,top.elev=7.04,width=10,roughness=1,pump=0,upstream=9,downstream=10,name='Allen to Gearheart')
	box[28,]=c(type=1,weir.elev=0,top.elev=6.95,width=10,roughness=1,pump=0,upstream=9,downstream=10,name='Allen to Gearheart')
	box[29,]=c(type=1,weir.elev=0,top.elev=6.67,width=10,roughness=1,pump=0,upstream=10,downstream=11,name='Gearheart to Hauser')
	box[30,]=c(type=1,weir.elev=0,top.elev=6.74,width=10,roughness=1,pump=0,upstream=10,downstream=11,name='Gearheart to Hauser')
	box[31,]=c(type=1,weir.elev=0,top.elev=7.02,width=10,roughness=1,pump=0,upstream=11,downstream=12,name='Hauser to Brackish')
	box[32,]=c(type=1,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=12,downstream=15,name='Brackish to Bay')
	box[33,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=2,downstream=7,name='Ox2 to Pilot')
	box[34,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=7,downstream=4,name='Pilot to T1')
	box[35,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=3,downstream=8,name='Ox3 to Cl')
	box[36,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=3,downstream=8,name='Ox3 to Cl')
	box[37,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=2,downstream=8,name='Ox2 to Cl')
	box[38,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=13,downstream=8,name='M-Box to Cl')
	box[39,]=c(type=2,weir.elev=0,top.elev=0,width=10,roughness=1,pump=0,upstream=11,downstream=8,name='Hauser to Cl')
